Geometry induced charge separation on a helicoidal ribbon 



Victor AtanascvQ 
Institute for Nuclear Research and Nuclear Energy, 
Bulgarian Academy of Sciences, 72 Tsarigradsko chaussee, 1784 Sofia, Bulgari^ 

Rossen Dandoloff 
Laboratoire de Physique Theorique et Modelisation , 
Universite de Cergy-Pontoise, F-95302 Cergy-Pontoise, Franc^ 

Avadh Saxena 
Theoretical Division and Center for Nonlinear Studies, 
Los Alamos National Laboratory, Los Alamos, NM 87545 USJw 

We present an exact calculation of the effective geometry-induced quantum potential for a particle 
confined on a helicoidal ribbon. This potential leads to the appearance of localized states at the 
rim of the helicoid. In this geometry the twist of the ribbon plays the role of an effective transverse 
electric field on the surface and thus this is reminiscent of the quantum Hall effect. 
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The interplay of geometry and topology is a recurring 
theme in physics, particularly when these effects mani- 
fest themselves in unusual electronic and magnetic prop- 
erties of materials. Specifically, helical ribbons provide 
a fertile playground for such effects. Both the helicoid 
(a minimal surface) and helical ribbons are ubiquitous 
in nature: they occur in biology, e.g. as beta-sheets 
in protein strucutres^l, macromolecules (such as DNA)0, 
and tilted chiral lipid bilayers^l Many structural motifs 
of biomolecules result from helical arrangements^: cel- 
lulose fibrils in cell walls of plants, chitin in arthropod 
cuticles, collagen protein in skeletal tissue. Condensed 
matter examples include screw dislocations in smectic 
A liquid crystaliP, certain ferroelectric liquid crystal^, 
and recently synthesized graphene ribbons. In particu- 
lar, graphene Mobius strips have been investigated for 
their unusual electronic and spin properties^. A helicoid 
to spiral ribbon transition^ and geometrically induced 
bifurcations from the helicoid to the catenoicP have also 
been studied. 

Graphene ribbons can be doped with charges. In this 
context, our goal is to answer the following questions: 
what kind of an effective quantum potential does a charge 
(or electron) experience on a helicoid or a helical ribbon 
due to its geometry (i.e., curvature and twist)? If the 
outer edge of the helicoid is charged, how is this poten- 
tial modified and if there are any bound states? Our 
main findings are: the twist lu will push the electrons 
in vanishing angular momentum state towards the inner 
edge of the ribbon and the electrons in non-vanishing 
angular momentum states to the outer edge thus creat- 
ing an inhomogeneous effective electric field between the 
inner and outer rims of the helicoidal ribbon. This is 
reminiscent of the quantum Hall effect; only here it is 
geometrically induced. We expect our results to lead to 
new experiments on graphene ribbons and other related 
twisted materials where the predicted effect can be ver- 



ified. In a related context we note that de Gennes had 
explained the buckling of a flat solid ribbon in terms of 
the ferroelectric polarization charges on the edges^. 
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FIG. 1: A helicoidal ribbon with inner radius £o and outer 
radius D. For £ = it becomes a helicoid. Vertical axis is 
along x and the transverse direction £ is across the ribbon. 

In order to answer the questions posed above, here 
we study the helicoidal surface to gain a broader un- 
derstanding of the interaction between quantum parti- 
cles and curvature and the resulting possible physical ef- 
fects. The properties of free electrons on this geometry 
have been considered beforeW The results of this pa- 
per are based on the Schrodinger equation for a confined 
quantum particle on a sub-manifold of K 3 . Following 
da Costa^ an effective potential appears in the two di- 
mensional Schrodinger equation which has the following 
form: 

Vcurv = — (M 2 — K) , (1) 

2m* 

where m* is the effective mass of the particle, h is the 
Plank's constant; M and K are the Mean and the Gaus- 
sian curvature, respectively. 
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To describe the geometry we consider a strip whose 
inner and outer edges follow a helix around the a;-axis 
(see Fig. [ljwith £ = 0). The surface represents a helicoid 
and is given by the following equation: 



f = x e x + £ [cos(tt>a;) e v + sm(ujx) e 2 



(2) 



where to — ^p 1 , L is the total length of the strip and n 
is the number of 27r-twists. Here (e x , e y ,e z ) is the usual 
orthonormal triad in K 3 and £ E [0,13], where D is the 
width of the strip. Let dr be the line element and the 
metric is encoded in 



u(%) 




-0,2- 



-0,4 



\dff = (1 + uj 2 ^)dx 2 + d£ 2 = hfdx 2 + h 2 2 d£ 2 , 



where hi = h\((,) = yl + w 2 £ 2 and h 2 = I are the Lame 
coefficients of the induced metric (from K 3 ) on the strip. 
Here is an appropriate place to add a comment on the 
helicoidal ribbon, that is a strip defined for £ S [£oj-D] 
(see Fig. [IJ. All the conclusions still hold true and all of 
the results can be translated using the change of variables 



t = Ho + s(D-Zo), 



se [0,1]. 



Here s is a dimensionless variable and one easily sees that 
for -> we again obtain the helicoid. 

The Hamiltonian for a quantum particle confined on 
the ribbon is given by: 



H 



1 



2m* hi 



d d 
dl hl dl 



d_J_d_ 
dx hi dx 



(3) 



Let us elaborate on the curvature-induced potential 
Vcurv ■ Since the helicoid is a minimal surface M vanishes 
and we are left with the following expression 



2m* 



-K 



2m* [l + w 2 C 2 ] : 



(4) 



Using Gauss' Theorema egregiunv^ the above potential 
can also be rewritten as 



2m* 



-K=- 



1 (d 2 hi 



2m* hi V <?£; 



(5) 



After rescaling the wave function ip i— > ^=^"0 (because 
we require the wave function to be normalized with re- 
spect to the area element dxdt;) we arrive at the following 
expression for the Hamiltonian: 



H 



2m* 



1 d 2 

h^dx 2 



+ V eff (H), (6) 



where the effective potential in the (transverse) £ direc- 
tion is given by: 



Veff(0 



2m* 



J_ f d 2 hi 
2hi \ «9£ 2 



IJL ( dhl 



(7) 



FIG. 2: The behavior of the potential U (£) for oj = 1 and 
ft 2 — 2m* — 1. Here the red curve corresponds to m = 0, the 
green curve to m = 1, and the yellow line to the approxima- 
tion given by Eq. (141. 



Note that in bent tubular waveguides^ and curved 
quantum strip waveguides^ the effective potential is lon- 
gitudinal. In the present case there is no longitudinal ef- 
fective potential. After insertion of hi 
effective potential becomes: 



Veff(0 



4m* (1 + w 2 £ 2 



1 



y/l + uj 2 C the 

u 2 e~ 



(8) 



This effective potential is of pure quantum-mechanical 
origin because it is proportional to h. Note that this 
expression is exact and is valid not just for small £: here 
no expansion in a small parameter has been used. 

Next, we write the time-independent Schrodinger 
equation as: 



d 2 



2m* <9£ 2 



Veff(0 



1> 



1 d 2 ^ 



2m* h\ dx 2 



Eip. (9) 



Using the ansatz: ip(x,t;) — 4>(x)f(£) we split the depen- 
dence on the variables and we get two differential equa- 
tions: 



h 2 d 2 (j){x) 



2m* dx 2 



E (f>(x), 



and 



where 



h 2 

2m* 



d 2 f(Q 
dt, 2 



u(0 = v eff (0 



Eq 
h 2 i(0' 



(10) 



(11) 



(12) 



With a solution <j)(x) = e " x of Eq. (10 1 we have 



h 2 

2m* ' 



where k x is the partial momentum in x-direction. Let us 
consider here the azimuthal angle around the x axis: lux 
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and the angular momentum along this axis: L x — — ^ • 
This operator has the same eigenfunctions L x <p(x) = 
hm<j>(x) as the operator in Eq. (10). The corresponding 
eigenvalues are Km. We conclude that the momentum k x 
is quantized 



tow, 



to € N. 



This is not surprising because of the periodicity of the 
wave function along x. Note that the value of the an- 
gular momentum quantum number determines the direc- 
tion the electron takes along the x axis either upward 
to > or downward to < 0. This situation is reversed for 
a helicoid with opposite chirality. 

Equation (11 1 represents the motion in the direction £ 



with a net potential 



u(0 



h 2 Ld 2 f 1 - 4to 2 



2m* 4 \(l+Lo 2 e) {l + uJ 2 £, 2 ) 2 



(13) 



which is depicted in Fig. [2j 

This potential is a sum of two contributions, an attrac- 



tive part: 



(i+^ 2 5 2 ) : 



and a variable part which is repulsive 



for to > 1 and attractive for to = (see Fig |2| . The ac- 
tion of this part for to ^ qualifies it as a centrifugal po- 
tential. It pushes a particle to the boundary of the strip. 
The finite size of the width D determines the cut-off of 
[/(£) and hence the probability of finding the particle is 
greatest near the rim of the helicoid. Since the behavior 
of the potential Z7(£) for a particle with to — qualifies 
it as a quantum anti-centifugal one, it concentrates the 
electrons around the central axis for a helicoid (or the 
inner rim for a helicoidal ribbon). Such anti-centrifugal 
quantum potentials have been considered before^. 

The behavior described above can be inferred using 
the uncertainty principle. Localized states must appear 
away from the central axis or the inner rim. Physically, 
one may understand the appearance of localized states 
away from the central axis using the following reason- 
ing: for greater £ a particle on the strip will avail more 
space along the corresponding helix and therefore the 
corresponding momentum and hence the energy will be 
smaller than for a particle closer to the central axis. 

We note that the separability of the quantum dynamics 
along x and £ directions with different potentials points 
to the existence of an effective mass anisotropy on the 
helicoidal surface. 

For the sake of simplicity let us approximate the po- 
tential Z7(£) given in Eq. (13 1 (for to = 1) by a straight 
line. The sole purpose of this approximation is to pin- 
point the basic distribution of the probability density. 
Assuming it to be linear (see Fig. [2| and starting from 
certain £o = a <C 1 



UaiO = (£>-£) 



Uo 



D 



U = */(£ = a). (14) 



The value of a can be determined from an area preserving 
condition ^Uq(D — a) — U{x)dx, where £o < D is 



the position from which we evolve the surface. When 
dealing with a helicoidal ribbon we must take ^ 0. 
After obtaining a result for this case we can easily obtain 
a result for the helicoid case by taking the limit £o ~- ► 0. 

Next we introduce a characteristic lengthscalc I in the 
problem 



r 



2m*\U \ 

h 2 {D~ a y 



x 
p 



2m* 



E 



DUo 
D-a 



where A is a dimensionless energy scale. After intro- 
ducing the dimensionless variable £ = — A — £// the 
Schrodinger equation for the radial part becomes 



d( 2 



C/K) = o, 



(15) 



with the following boundary conditions: /(—A — £o/Z) = 
/(—A — D/l) = 0. This form of the equation is valid for 
Uo > as is the case for m^O. 

For to = we have a negative Uo = —\Uo\ which 
requires the introduction of the dimensionless variable 
£ = — A + £/Z and the corresponding equation is given 



by (15 1, only in this case the boundary conditions are 
/(-A + £ /0 = /(-A + J D/0 = 0. 

Let us assume that the ratio D/l 3> 1 then the solu- 
tions, i.e. the wave functions, of Eq. (151 coincide with 
the Airy function, that is /(£) = const Ai(£), and the 
boundary condition /(—A ± £o//) = (the upper sign 
corresponds to m — and the lower to to ^ states) 
gives the quantized energies 



E n (m) = U {m) 



D 



D-a 



A„± 



£o 



I J 2m*l 2 ' 



where A„ are the zeroes of the Airy function Ai(— A„) = 
0. Let us list the first three of them: (Ai,A2,A 3 ) = 
(2.338,4.088,5.521). Here we have taken account of the 
case when the interior of the helicoid is cut at a distance 
£o from the axis, that is the ribbon case. The helicoid 
case is obtained after setting £ — > 0. 

For the vanishing angular momentum state we have 
Z7o(0) < and the energy spectrum starts at a negative 
value (Fig. [2J, that is we have a bound state. The prob- 
ability amplitude has a node at £o in the ribbon case or 
at the origin for the helicoid case. The evolution along £ 
starts at the corresponding zero of the Airy function and 
evolves in the positive direction where the Airy function 
vanishes. For non- vanishing angular momentum states 
we have Uq (to) > and the energy spectrum is positively 
valued (Fig. [2]). The evolution of the corresponding so- 
lutions along £ starts at the corresponding zero of the 
Airy function and evolves in the negative direction where 
the Airy function is oscillatory as one would expect for 
a confined positive energy spectrum. The observation 
that the to ^ states at £o = X n l have the same energy 
E n (m) — Uo(m)D/(D — a) for all n leads us to believe 
that this is a particular positive energy oscillatory state 
whose wavelength fits D(l - £ /Z) « D(l/D > 1). 
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We would like to conclude with the observation that 
the electric dipole moment for the (to = 0,Ai) bound 
state (also the ground state for this geometrical config- 
uration) is non-zero due to the anisotropic distribution 
of the probability density along £. Indeed, suppose we 
consider a ribbon doped with a uniform surface charge 
density a, then the electric dipole vector p = p x e x + p^e^ 
in the moving coordinate system [e x ,e^,es = e x x eg) 
will have non- vanishing x and £ components: 

Px = , Pe = — crl fj n , (16) 

where the total charge is 

/■27t/oj pD 

Q= dx' <r\i/)(a/,?)\ 2 d? 

JO J£ 

and 




Here the upper sign corresponds to to = and the 
lower to to 7^ states. For £ /l = 0.1 and D/l = 10 we 
summarize the values of (3 n in the following table 





n= 1 


n = 2 


n — 3 


n = 10 


m = 


0.816 


1.822 


2.829 


3.605 


m ± 


2.712 


2.451 


2.299 


1.783 
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